Mixing of Quantum Walk on Circulant Bunkbeds* 
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Abstract 

We give new observations on the mixing dynamics of a continuous-time quantum 
walk on circulants and their bunkbed extensions. These bunkbeds are defined through 
two standard graph operators: the join G + H and the Cartesian product G ® H of 
graphs G and H. Our results include the following: 

• The quantum walk is average uniform mixing on circulants with bounded eigen- 
value multiplicity. This extends a known fact about the cycles C n . 

• Explicit analysis of the probability distribution of the quantum walk on the join 
of circulants. This explains why complete partite graphs are not average uniform 
mixing, using the fact K n = K\ + K n -\ and iCi,...,n = K n + . . . + K n . 

• The quantum walk on the Cartesian product of a m-vertex path P m and a circulant 
G, namely, P m © G, is average uniform mixing if G is. This highlights a difference 
between circulants and the hypercubes Q n = P2 ffi Q n -i- 

Our proofs employ purely elementary arguments based on the spectra of the graphs. 
Keywords: Quantum walks, Circulant graphs, Average mixing, Join, Cartesian product. 

1 Introduction 

The study of continuous-time quantum walks on graphs has important potential applications 
in quantum computation |T£]. First, as an algorithmic technique, it was used to devise 
efficient quantum search algorithms with considerable speedup over classical algorithms [7]. 
Second, it may provide a simpler physical implementation of a quantum computer, given 
that there is an abundance of physical processes that simulate quantum walk on graphs 
[S]. In the physics literature, continuous-time quantum walks is mainly studied over infinite 
constant-dimensional lattices, such as the one-dimensional line (see [HJj, Chapters 13,16). 
On the other hand, the study of random walks on general graphs is a topic of broad interest 
in the mathematics and computer science community (HI El • 
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In this paper, we study the mixing dynamics of continuous-time quantum walks on circu- 
lant graphs. More particularly, we consider the average or limiting probability distribution 
of a quantum walk. This notion was introduced in and is the quantum analogue of a 
stationary distribution of classical random walks. On the circulant graphs, our goal was 
to characterize the graphs for which the continuous-time quantum walk reaches (almost) 
uniform average probability distribution. It was previously known that cycles are near uni- 
form mixing, whereas the complete graphs and hypercubes are not (see (HE])- Our other 
goal in this paper is to discover graph theoretic structures that may explain this polarized 
phenomena. 

First, we show that circulants with bounded eigenvalue multiplicity are almost uniform 
mixing. This generalization explains why cycles are uniform mixing. Second, we consider 
bunkbed graphs constructed using the join and the Cartesian product operators. By ana- 
lyzing the join of two circulants, we observe an interesting mixing phenomena on the cone 
K\ + G of a circulant G, that is dependent on the density of G. If the quantum walk starts 
on K\, a dense graph G repels the probability away from the copy of G. This explains why 
the limiting distribution of a quantum walk on the complete graph is not near the uniform 
distribution. We extend this investigation to the homogeneous join of circulants, namely, 
G + . . . + G, for a circulant G. We show that this bunkbed graph is uniform mixing if G is 
uniform mixing and the join is over a constant number of copies of G. A corollary of this 
transference property explains the non-uniform mixing of the complete multipartite graphs 
K n + ... + K n . 

We also analyze bunkbed graphs obtained from the Cartesian product P m © G of a 
path P m and a circulant G. On this bunkbed structure, we observe another transference 
property: the quantum walk on P m © G is uniform mixing if it is uniform mixing on G and 
the path is of constant size. This highlights a striking difference with the hypercube Q n , 
since the hypercube is also a bunkbed Q n = P2 © Q n -i, but it is known that they are not 
uniform mixing [Tl]. It is interesting to note that both classes of graphs are group-theoretic 
circulants (see (Hj), since our circulants are the Z n -circulants while hypercubes are the 
(Z2) n -circulants. This suggests a group theoretic investigation into the mixing phenomena 
of generalized circulants, which we leave for future work. 

In this paper, we focus exclusively on continuous-time quantum walks. We refer the 
reader to J2j for a survey of other models of quantum walks. As a final remark, we 
mention that most of the graphs we consider have the standard stationary distributions in 
the classical random walks where the limiting probability of a vertex is proportional to its 
degree (!]. 

2 Preliminaries 

We consider simple, undirected graphs that are connected, and mostly regular. For a graph 
G = (V,E), let Aq be the adjacency matrix of G, where Aq[j, k] = {(j,k) £ Ej. Here and 
throughout, we will use to denote the characteristic function of a logical statement 
that is, 1 if ^ is true, and if it is false. The set of eigenvalues of Aq is denoted Sp(G), and 
the (algebraic) multiplicity of an eigenvalue A is denoted m(A). The spectral type r(G) of 
a graph G is the number of distinct eigenvalues of the adjacency matrix Aq of G. We will 
denote the maximum (algebraic) multiplicity of any eigenvalue of graph G by n(G). Some 
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Figure 1: Examples of circulants of order 8. From left to right: (i) the empty graph K%. 
(ii) the cycle C%. (iii) a strongly-regular circulant: clique minus a perfect matching, (iv) 
the complete graph Kg. 

of the families of graphs that we will consider include the complete multipartite graphs 
where there are m partitions with a partition size of n, the cycles C n and paths P n , 
and the hypercubes Q n . Relevant background on graphs and their spectral properties can 
be found in jSJ. 

A graph G is called circulant if its adjacency matrix Aq is circulant. A circulant matrix 
A is specified by its first row, say (do, a\, . . . , a n _i), and is defined as Aj^ = a k _j ( mod n ), 
where j,k £ 7L n . Here 7L n denotes the group of integers {0, ... ,n — 1} under addition 
modulo n. Note that ao = 0, since our graphs are simple, and cij = a n -j, since our graphs 
are undirected. Connectivity is guaranteed if the greatest common divisor of n and all 
indices k, for which a k = 1, is one. Alternatively, a circulant graph G = (V,E) can be 
specified by a subset S C Z n , where (j, k) £ E k — j £ S. In this case, we write G = {S). 
We will assume that S is closed under taking inverses, namely, if d & S, then — d £ S. 
Figure H contains some examples of circulant graphs. 

It is known that circulant graphs G are diagonalizable by the Fourier matrix F defined 
as Fj k = n -1 / 2 ti^ fc , where ui n = exp(27ri/n). In fact, the eigenvalues of Aq are 

n-l L(»-1)/2J /„ .,v 

A J = ^a fc ^ fc = 2cos(^-) + [n even] a n/2 (-r/. (1) 

k=l k=l \ n J 

A continuous-time quantum walk on a graph G = (V, E) is defined using the Schrodinger 
equation with the real symmetric matrix Aq as the Hamiltonian (see [7]). If \ip(t)) £ C^l is 
a time-dependent amplitude vector on the vertices of G, then the evolution of the quantum 
walk is given by 

m)) = e- itAG \m), (2) 

where i = and \ip(0)) is the initial amplitude vector. We usually assume that |^(0)) is 
a unit vector, with (x\ijj(0)} = \x = 0], for some vertex 0. The amplitude of the quantum 
walk on vertex j at time t is given by ipj(t) = 0I^(^)); while the probability of being on 
vertex j at time t is Pj(t) = |^j(i)| 2 . The average (or limiting) probability of being on 
vertex j is defined as 

1 f T 

Pj = Jim t= / Pj{t) dt. (3) 

This notion appeared in in the context of discrete-time quantum walks. The limiting 
probability distribution of the quantum walk will be denoted P. 



3 




Figure 2: Some circulants with bounded eigenvalue multiplicity. From left to right: (i) the 
wheel a cycle plus a perfect matching, (ii) a double-loop circulant (±1, ±2) of order 7. 



Definition 1 (Average Uniform Mixing) 

The average mixing of a continuous-time quantum walk on a graph G = (V, E) is called 
uniform ifpj = 0(1/\V\), for each vertex j of G. 

Remark Note that in the above definition, we only require that each limiting probability be 
linearly proportional to the uniform probability value. This is less stringent than requiring 
that the quantum walk achieves exactly uniform probability distribution (see HI ) ■ When 
the graph G is not regular, the limiting probability distribution P may depend on the 
initial state \ip(0)). We will specify carefully the effect of the initial states in these cases, 
but suppress this dependence for vertex-transitive graphs. 



3 Mixing and Bounded Multiplicities 

Theorem 1 Let G be a circulant graph. If fJ.(G) is bounded, then the continuous-time 
quantum walk on G is average uniform mixing. 

Proof Let n be the order of G and let A be the adjacency matrix of G. Since |0) = 
E n k Zl^\F k ), if |V(0)) = |0), we have ty(t)) = e~ aA \m) = ^Ek=J e " <tt W This 
yields (j\i;(t)) = ~ Zk=o e- itx ^ k . Thus, 

P f t ) = J_ y e -ft(A k -^-(*-o = I + J_ y e -»t(A fc -A,)^ M . (4 ) 

n l z — ' n n 2 ^-^ 

k,e k^t 

Using the above, the average (limiting) probabilities are 



Pj 

n 



<- i E (t») <- m 

xes P (A) v 7 v 7 

So, if n{G) = 0(1), we have uniform mixing. □ 

Theorem 2 The continuous-time quantum walk on a constant- degree n-vertex circulant of 
the form (1, n/ki, . . . , n/k^), where, for each j = 1, . . . , d, kj is a constant which divides n, 
is average uniform mixing. 

Proof The eigenvalues of G are given by 
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for j = 1, . . . ,n — 1. Since the sum Yle=i cos (2tu/^£) can have at most Yle=i = 0(1) 
distinct values, each eigenvalue must have a constant multiplicity. By Theorem ^ we have 
the claimed result. □ 

Corollary 3 The continuous-time quantum walk on the 3-regular circulant "wheel" V n = 
(l,n/2) of even order n is uniform mixing. 

4 Mixing on Join Bunkbeds 

In this section, we study a circulant bunkbed structure obtained by the join of circulants. 
Formally, the join G + H of two graphs G and H is defined as to satisfy G + H = G U H 
(see ^E])- It is easy to see that this is a graph obtained by connecting each vertex of G to 
each vertex of H, while maintaining the internal structures of G and H. For a graph G, 
the cone of G will denote the graph K\ + G. 

Lemma 4 Let G and H be circulants of degrees k and t, respectively. Suppose that the 
eigenvalues of G and H are k = Ho > fj,% > . . . > fA\G\-i an d & = u o > v\ > • • • > v \H\-ii 
respectively. Then, the eigenvalues and (orthonormal) eigenvectors of G + H are found in 
three separate sets {(fi a ,\z^)) : 1 < a < |G| — 1}, {{vb,\ z b)) '■ 1 — & — 1-^1 ~~ an ^ 
{(^± ; |<z=fc))}j where, for x = 0, . . . , |G| \H\ — 1, we have 

(x\ Z C) = -4^[xeG], a=l,...,\G\-l (7) 



\G\ 

(x\z?) = u\ H{ {x e Hj, b = l,...,\H\-l (8) 



(x\z±) = (/3±)^, (9) 
■L>± 



where (3± = (X± — k)/\H\, L± = y\G\ + \H\j3±, and X± are the roots of X 2 - (k + £)X 
(\G\\H\ - kl) = 0. 



Proof Note that the adjacency matrix of G + H is given by 



A 



Ag J\G\x\H\ 
J\H\x\G\ Ab 



(10) 



It is easy to see that \z^} are eigenvectors of A with eigenvalues fx a , for a = 1, . . . , |G| — 1, 
and \z?) are eigenvectors of A with eigenvalues Uf,, for 6 = 1,..., \H\ — 1. The last two eigen- 
vectors are obtained by noting that the eigenvectors have the form [a ... a b ... b\ T . 
This gives the equations ka + b\H\ = Xa and lb + a\G\ = Xb, whose solutions yield the 
eigenvalues X± = \{{k + I) 2 ± \/A), where A = (A; — £) 2 + 4|G||i7|, and eigenvectors with 
a = land&=(A±-fc)/|.ff|. □ 

Theorem 5 Suppose that G and H are circulants of degrees k and I, respectively. Let 
A = (k — €) 2 + A\G\ \H\ and X± = l/2[(fc + £) 2 ± y/A]. Consider a continuous-time quantum 
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walk on G + H starting at some vertex of G. Let p x {G) denote the limiting probability of 
x £ G over the subgraph G. Assume that 

\-?(Sp(G)\{k})u(Sp(H)\{£}). (11) 

Then, the limiting probabilities of the vertices of G + H are 

pAG + H) = !(p x {G)-^+^f^-^^hxeGj + ^{xeHj (12) 

Proof Let the initial state be \ip(0)) = |0) where the quantum walk starts at a vertex of G. 
By Lemma 0] we have 

1 |GM 1 



and, thus, 



1 ^ p~^± 



\4>V» = -m ^ ^l+EV- |z±) - (14) 

VI G I =1 ± * 



The amplitude on vertex a; at time i is given by 



-i IGI- 1 -itx ± 

ww> = fi E e "^-ici + E V-* 631 ' (15) 

where /3± = (A± — and we obtain 

1 f T 1 f 6 lxeH} \ 2 

p x =lun-l dt \(xm))\ 2 =P x (G)- |^ + E(^^I » 



where p x {G) is the limiting probability on the subgraph G. After some calculations, we get 

= m(m _ V)' 5 (ft) =z ' (17) 

which completes the stated claim. □ 

The single theorem above implies the following various known and new facts about mixing 
on the family of complete and related graphs. First, we obtain a perfect uniform mixing 
behavior on Ki, but not on K n , for n > 2. 

Corollary 6 The continuous-time quantum walk on K2 is average exactly uniform mixing. 

Proof By Theorem 03 we have \G\ = \H\ = 1 and k = I = 0. Thus, A = 4, and therefore, 
p =p 1 = 1/2. □ 

Corollary 7 [2] The continuous-time quantum walk on K n is not average uniform mixing, 
as n — > 00. 
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Proof By Theorem^ we have K n = K\ + K n -\. We have \G\ = 1, \H\ = n — 1, fc = 0, and 
£ = n-1. Then, A = (n-l) 2 +4n, withp = l-2n/A ~ 1 andpj = 2/ A ~ 0, as n -> oo. □ 

Next, we consider the cone of circulants. The following corollary provides a simple expla- 
nation why K n is not average uniform mixing, for large n; it is because K n is a cone of a 
dense circulant. 

Corollary 8 The continuous-time quantum walk on the cone of any circulant C , namely, 
K\ + C, is not average uniform mixing. 

Proof Let C be a ^-regular circulant of order n. By Theorem|21 we have \G\ = 1, \H\ = n, 
k = 0. Then, p = 1 - (1/2)[1/(1 + (^/2) 2 /n)]. Thus, p = 0(1), regardless of £. D 

Homogeneous Joins of Circulants Consider the unbounded m-fold homogeneous join 
of a circulant G, namely, G^ +m ^ = G + . . . + G, where there are m terms in the summation. 
The following theorem shows that the uniform mixing property of G transfers into its 
unbounded homogeneous join if m is a constant. 

Theorem 9 Let G be a circulant of order n. Let m > 2 is a constant and n > 2Ao(G). In 
the continuous-time quantum walk, G^~ m ' = EaLi G is average uniform mixing if G is. 

Proof The adjacency matrix of G^ +m ^ = E&Li G * s given by 

A = I m <g> G + K m ® J n , (18) 

where I m is the mxm identity matrix, K m is a complete graph on m vertices, and J n is 
the n x n all-one matrix. Since G is a circulant, both summands share the same set of the 
following orthonormal eigenvectors 

{\F jik ) = |Fj m) >® |i> } > : 0<j<ro-l,0<fc<n-l}, (19) 

where | F^ ) denotes the j-th column of the mxm Fourier matrix, and similary for | F^ ) . 
Let \k(G), for < k < n— 1, be the eigenvalues of G in descending order. The corresponding 
eigenvalues of C?( +m ) are given by 

r A (G) + (m-l)n ifj = jfe = 
X jt k = { A (G) - n if j / and k = (20) 

I A fc (G) if 3,k^0 

If |V(0)> = |0) ® |0) then |^(0)) = l/v^E^E*^ l*i,k>- Thus > 

^ m n— 1 

|^)) = ^^^e-^|F,, fe ). (21) 
v j=0 fc=0 

Thus, for x G Z m and y 6 Z n , we have 

^ m n— 1 

^(i) = (x, = V V e-^.» exp 
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Note that the three types of eigenvalues of G( +m ) are mutually distinct, since 

A (G)-n < A fc (G) < Ao(G) + (m - l)n. (23) 

Therefore, we have 



1 

Px,y 



mn 



1 / m — 1 

< 



(mn) 2 \ 2 



V 2 ;j 



(^') (2 1.) 

mn • 



since m is a constant. □ 

The above theorem also explains why the complete graph Kn is not uniform mixing, since 
Kn can be viewed as a homogeneous m-fold join of K N / m , for some constant m that divides 
N. The theorem also implies the following claim about the multipartite complete graphs. 

Corollary 10 The continuous-time quantum walk on the complete multipartite graph 
is not average uniform mixing if m > 2 is a constant. 

Proof Since a continuous-time quantum walk is not average uniform mixing on the empty 
graph K n and = K n + . . . + K n , we have our claim. □ 

5 Mixing on Cartesian Bunkbeds 

In this section, we consider a circulant bunkbed structure obtained by the Cartesian product 
i-2 © C, where C is a circulant graph. 



Lemma 11 Let G be a circulant of degree d and order n, whose eigenvalues are d = uq > 



Hi > ... > A„_i. Then, the eigenvalues of P2 © G are = jij ± 1 with the following 



(orthonormal) set of eigenvectors 

\zf) = \±)®\ Zj ), (25) 

where \±) = ^(|0) ± |1)) and (x\zj) = (1/ y / jG\)uj jx , for j,x G [n], with oj = e 2ni / n . 

Proof Note that the adjacency matrix of P2 © G is given by P2 © I n + I2 © Ag- Since P2 
is a circulant, both P2 © I n and I2 © Aq are simultaneously diagonalizable by \z^). This 
implies the stated claim on the spectra of P2 © G. □ 

Theorem 12 Let G be a circulant of order n. In the continuous-time quantum walk, P2©G 
is average uniform mixing if G is. 

Proof Assume that |V>(0)) = |0) © |0). Thus, |V(0)> = £± ^|±) © J2j and 

\ m ) = e~ itA \m) = 4= E e-*>?\±) © |^->. (26) 
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Figure 3: Examples of circulant bunkbeds. (a) The double bunkbed: i-^ffiCs- (b) The Join 
bunkbed: C 8 + C 8 . (c) The Cone: K x + C 8 . (d) The Cartesian bunkbed: P 4 © C 8 . 



This implies that 

<6,x^(t)> = -^=^e-^ ±1 )(6|±)( a; |z j ) (27) 

= —^2e- it ^ ±1 \±l) b uj jx (28) 

= -^e- ia ^^^e- i * (±1) (±l) 6 (29) 
n j ± 

I ^ e -^^ :c [(l - 6) cos(t) + b(-i sin(t))] (30) 



n 



Letp M (t) = |(6,x|^(t))| 2 . Thus, 

P6,*(*) = ^ S^ 1 " 6 ) cos2 ^) +bsm 2 (t)]e- it ^- x ^LO ( - j -V x . (31) 



Note that po, x (t) + pi, x (t) = p x (t), where p x (t) is the (instantaneous) probability on vertex 
x at time t of a quantum walk on G alone. Then, 

Po,x 4E w(, ^ T lim J f T cos\t)e-^-^ At = \ p x , (32) 
j,k 70 

since limT_ >00 y Jq cos 2 (t)e~ ltA dt = \ [A = 0]. Similarly, we obtain = ^p^.. This 

yields the claim. □ 

Corollary 13 XTje continuous-time quantum walk on a Cartesian bunkbed P2 © G, where 
G is a d-degree circulant of the form (l,n/fei, . . . ,n/kd-i), where d and k±,..., kd-i are 
constants, is average uniform mixing. 
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Circulant Cylinders To extend our Cartesian bunkbeds over paths with more than two 
vertices, we provide, for completeness, an analysis of the quantum walk on paths. This 
problem is well-known in the physics literature, but is normally done on the infinite paths 
using different techniques The eigenvalues Xj and eigenvectors \Qj) of the path P m 

(see ^7|), for j = 1, ... , m, are defined as 

x ' = 2 ™Gstt) (33) 

1 / J X7T \ 

(x\Qj) = — = sin , x = 1, .... m (34) 

V y/(m + l)/2 \m+l) y ' 

If the quantum walk starts with the initial state 1^(0)) = where the basis states are 
|1), . . . , \ m), then Thus, we have 

m —it\ / ■ \ 



- ^(m+l)/2 



Since P m has m distinct eigenvalues, the limiting probabilities are given by 

4 ST ■ 2 ( 3^ \ . 2 ( jx-K \ 

Px = 7 tttt > sin sin . (36) 

x (m+l) 2 p^ \m + l) \m+l) K ' 

Note that, since sin 2 (t)di = vr/2, we get an upper bound of 

4 • 2 / J 71 " 



Px 



< 



(m + l) 2 ^— ' V + 1 

v ; i=i 



< 7 ^— f / sin 2 (t) dt + , - - ) (38) 

" (m + l)n\J (m + l)J 

^ T^Tu + ( lu2 =°(-)> ( 39 ) 
(m + 1) (m + l) z \m/ 

which implies that the quantum walk on P m is average uniform mixing. 

The eigenvalues of a circulant cylinder T = P m © G, where G is a circulant of order n, 
are given by 

Xj t k = (J>j + ^fc, where 1 < j < m, < k < n — 1, (40) 

where //j = 2 cos(j7r/(m + 1)) and are the eigenvalues of P m and G, respectively. Since 
the adjacency matrix of T is defined as P m ® I n + Im ® G, the eigenvectors of T are 

\T j>k ) = \Qj) <g) \F k ), where 1 < j < m, 0<k<n-l, (41) 

where \Qj) and ji*}} are the eigenvectors of P m and the circulant G, respectively. Recall 

for 



that (x\Qj) = for 1 < j,x < m, and (y\F k ) = ^exp^J, 

< k,y < n — 1. 
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If the initial state is \ip(0)) = |1) <S> |0), we have 

m n—1 

1^(0)) = £<WlQi> ® ^(F fc |0)|F fc ). (42) 

j=l fc=0 

The adjacency matrix of PmSG is given by A = P m ® I n + I m ®G ', where the two summands 
commute with each other. Thus, e~ ltA = e ~ tt ( Pm ® In ') e~ lt ( Im ® GS> , and 

m n—1 

\m) = ^(Q^'^W ® E^|0)e-^|F,). (43) 

j=l k=0 

The amplitudes of \tp(t)) at vertex x on the path P m and vertex y within the circulant G is 
given by 

m n—1 

(x,y\m) = E^^I^^^I^E^^I^X^IO) (44) 
i=i fc=o 

Corollary 14 Let G be a Cartesian product P m ®C , where C is a circulant of order n. The 
continuous-time quantum walk on G is uniform mixing if m is constant or n is constant. 

6 Conclusions 

It was known that a continuous-time quantum walk is uniform average mixing on the cycles 
C n , but is not uniform average mixing on the complete graphs K n and on the hypercubes 
Q n . Our goal in this work was to provide a graph-theoretic explanation for this polarized 
phenomena. 

First, we extend the phenomenon of the cycles, by showing that uniform mixing is 
achieved on circulants with bounded eigenvalue multiplicity. We also gave other explicit ex- 
amples of circulants meeting this criteria. Second, we consider two graph-theoretic bunkbed 
structures over circulants in order to study the non-uniform mixing on K n and Q n . Our 
analysis on the join bunkbed sheds some light on the non-uniform mixing of the complete 
multipartite graphs (which includes K n ). Our analysis of the Cartesian bunkbed of circu- 
lants highlights a difference between the Z„-circulants and the (Z2) n -circulants (see (Hj). 
We leave a similar investigation of general group-theoretic circulants and Cayley graphs for 
future work (see [TT]). 
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